DISPERSIVE ESTIMATES USING SCATTERING THEORY FOR 
MATRIX HAMILTONIAN EQUATIONS 



(N 
3 



JEREMY L. MARZUOLA 

ON 

Abstract. We develop the techniques of [KriSchl and [ESI in order to derive dispersive 
estimates for a matrix Hamiltonian equation defined by linearizing about a minimal mass 
soliton solution of a saturated, focussing nonlinear Schrodinger equation 

iu t + Au + P(\u\ 2 )u = 
u(0,x) = u (x), 

in R 3 . These results have been seen before, though we present a new approach using 
scattering theory techniques. In further works, we will numerically and analytically study 
p j the existence of a minimal mass soliton, as well as the spectral assumptions made in the 

analysis presented here. 

•5 

, ^, 1. Introduction 

In this result, we develop the dipsersive estimates used to prove stability of solitons for 
^ | a focussing, saturated nonlinear Schrodinger equation (NLS) in R x WL d : 

J£j iut + Au + P(\u\ 2 )u = 

u(0,x) = u (x), 

where j3 : R — > R, /3(s) > for all s e R, /3 has a specific structure outlined in the following 

ON definitions: 
O 

j>! Definition 1.1. Saturated nonlinearities of type 1 are of the form 

(1.1) /*(*) = ^T^ET, 

03 1 + s 2 

where p > 2 + | and 4 > g > /or d > 3 and oo>p>2+|>|>g>0 /or d < 3. 
Definition 1.2. Saturated nonlinearities of type 2 are of the form 



(1.2) /?(*) 



where ^ > q > 0, d > 2. 



2-£ J 
S 2 



Remark 1.1. In k>£/i cases, for \u\ large, the b ehav ior is L 2 subcritical and for \u\ small, 
the behavior is L 2 supercritical. For Definition p is chosen much larger than the L 2 

l 
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critical exponent, | in order to allow sufficient regularity when linearizing the equation 
about the soliton. 

In the sequel, we assume that Uq € H l and \x\uq G L 2 , or in other words, Uq has finite 
variance. For initial data with this regularity, from the spatial and phase invariance of 
NLS, we have many the following conserved quantities: 

Conservation of Mass (or Charge): 



Q(u) = - j \u\ 2 dx = - j \u \ 2 dx, 



and 

Conservation of Energy: 

E{u)= [ \Vu\ 2 dx- [ G{\u\ 2 )dx= [ \Vu \ 2 dx- [ G(\u \ 2 )dx, 

JR d JR d JR d JR d 

where 

G(t) = [ (3(s)ds. 
Jo 

We also have the pseudoconformal conservation law: 
(1.3) \\{x + 2itV)u\\ 2 L2 -4t 2 [ G{\u\ 2 )dx = \\x<j)\\ 2 L 2 - [ 6{s)ds, 

JR d JO 

where 

9( s )= [ (4(d + 2)G(|w| 2 ) - 4df3(\u\ 2 )\u\ 2 )dx. 

J~R d 

Note that (x + 2itV) is the Hamilton flow of the linear Schrodinger equation, so the above 
identity shows how the solution to the nonlinear equation is effected by the linear flow. 

Detailed proofs of these conservation laws can be arrived at easily using energy estimates 
or Noether's Theorem, which relates conservation laws to symmetries of an equation. Global 
well-posedness in L 2 of (NLS) with f3 of type 1 or 2 for finite variance initial data follows 
from standard theory for L 2 subcritical monomial nonlinearities. Proofs of the above results 
can be found in numerous excellent references for (NLS), including |Caz] and [SulSulj . 
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2. Soliton Solutions 

A soliton solution is of the form 

u(t,x) = e iXt R x (x) 

where A > and R\(x) is a positive, radially symmetric, exponentially decaying solution 
of the equation: 

(2.1) AR X -XR X + (3(R X )R X = 0. 

With nonlinearities of type 1 or 2, soliton solutions exist and are known to be unique. 



Existence of solitary waves for nonlinearities of the type presented in Definitions 1.1 and 



1.2 is proved by in |BerLion] by minimizing the functional 
with respect to the functional 

V{u) = J [G(\u\ 2 ) - ^\u\ 2 ]dx. 

Then, using a minimizing sequence and Schwarz symmetrization, one sees the existence 
of the nonnegative, spherically symmetric, decreasing soliton solution. For uniqueness, 
see |McCleod] . where a shooting method is implemented to show that the desired soliton 
behavior only occurs for one particular initial value. 

An important fact is that Q x = Q(R X ) and E x = E(R X ) are different iable with respect 
to A. This fact can be determined from the early works of Shatah, namely [Shatahl] . 



[Shatah2] . By differentiating Equation (2.1), Q and E with respect to A, we have 



d x E x = —Xd x Q x . 

Numerics show that if we plot Q x with respect to A, we get a curve that goes to oo 
as A — > 0, oo and has a global minimum at some A = Ao > 0, see Figure [Tj We will 
explore this in detail in a subsequent numerical work |Mar-num] . Variational techniques 
developed in |GrilShaS~tr] and |ShatStri] tell us that when 5(A) = E x + XQ X is convex, 
or 5"(X) > 0, we are guaranteed stability under small perturbations, while for 5"(X) < 
we are guaranteed that the soliton is unstable under small perturbations. We will explore 
the nature of this stability in another subsequent work, |Mar-nonlin] , where we study the 
full nonlinear problem. For a brief reference on this subject, see |SulSul] . Chapter 4. For 
nonlinear instability at a minimum, see |ComPel] . For notational purposes, we refer to a 
minimal mass soliton as R m in- 
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Figure 1. Plots of the soliton curves (Q(X) with respect to A) for a sub- 
critical nonlinearity (d — 1, p — 3), supercritical nonlinearity (d — 3, p — 3), 
critical nonlinearity (d = 1, p = 5), saturated nonlinearity of type 1 ( p — 7, 
q = 3) in E, saturated nonlinearity of type 1 in 3d (p = 4, q — 2), saturated 
nonlinearity of type 2 in IR 3 (g = 2). The curves for the monomial nonlinear- 
ities are found analytically, while the curves for the saturated nonlinearities 
are found numerically. 
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3. Linearization about a Soliton 

Let us write down the form of (NLS) linearized about a soliton solution. First of all, we 
assume we have a solution ip = e lXt (R\ + <fi(x,t)). For simplicity, set R = R\. Inserting 
this into the equation we know that since R is a soliton solution we have 

(3.1) i(0) t + A(0) = -f3(R 2 )<p-2P\R 2 )R 2 Re(<P) + 0(<p 2 ), 

by splitting up into its real and imaginary parts then doing a Taylor Expansion. Hence, 
if (p = u + iv, we get 

(3.2) d t ( I ) = n 1 u 

where 

(3.3) H 
where 



v \ v 



L_ 
-L+ 



and 



-A + A - p{R> 



L+ = —A + A — f3(R x ) - 2P'(R 2 X )R 2 X . 



Definition 3.1. A Hamiltonian, TC, is called admissible if the following conditions hold: 

1) There are no embedded eigenvalues in the essential spectrum, 

2) The only real eigenvalue in [—A, A] is 7 

3) The values ±A are not resonances. 



Definition 3.2. Let (NLS) be taken with nonlinearity (3. We call (3 admissible if there 
exists a minimal mass soliton, R m in, for (NLS) and the Hamiltonian, 7i, resulting from 



linearization about R m i n is admissible in terms of Definition 3. 1 



The spectral properties we need for the linearized Hamiltonian equation in order to prove 



stability results are precisely those from Definition 3^ Notationally, we refer to Pd and 
P c as the projections onto the discrete spectrum of Tt and onto the continuous spectrum 
of 7i respectively. Analysis of these spectral conditions will be done both numerically and 
analytically in |Mar-spec] . 



4. Main Results 

We derive the existence and important properties of distorted Fourier bases 0g of non- 
self-adjoint matrix Hamiltonians, and hence a distorted Fourier transform, for a general 
class of matrix Hamiltonians. Let S be the Schwartz class of functions. Then, we have the 
following results: 
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Theorem 1. Given an admissible Hamiltonian 7i, and the projection on the continuous 
spectrum of TC, P c , for initial data 6 S, we have 

\\j tH Pc<j>\\L°° <H. 

Let us define the space 

L l)M = {fe ^IHO^/COIIli < oo, N = 0, 1, . . . ,2M}, 
with norm || ■ \\li,m defined in the standard fashion. 

Theorem 2. Let TC be an admissible Hamiltonian as defined above. Assume ip £ L 1, and 

(4.1) ^9^(0) = 0, 

for multi-indices a, (3 such that \a\ + \(3\ = 0, 1, 2, ... , 2M , where 



Then, 

(4.2) \\e- c ^e itn Pj\\ L ^ < Cr$- M m \ L i,M, 

for any c > 0. 

It should be noted that similar estimates were proven in the works jESlj and [BouWaJ, 
where in the first the techniques used were more along the lines of resolvent estimates and 
in the second the fact that the nonlinearities of interest were of even integer powers was 
crucial to the argument. Here, we take an approach similar to that of scattering theory 
as presented in |Ho2j . Scattering theory is related to a resolvent approach most certainly, 
though there are certain benefits to the method we thought would be of general interest. 
Note, these dispersive estimates are essential for the forthcoming argument in |Mar-nonlin] , 
where perturbations of minimal mass solitons are analyzed. 

5. General Distorted Fourier Basis Theory 



We present here a review of combined results from |Agmon and |Ho2j . Chapter 14. Both 
presentations are valid for operators of the form 

(P(D) + V(x,D))u = 0, 

where P{D) is a self-adjoint, constant coefficient differential operator and V(x,D) is a 
short range, symmetric differential operator. The perturbation V(x,D) is defined to be 
short range in order to say that 

lim R(z) = R^iz) 

exists in the uniform operator topology of B(L 2,S , 7i 2 ,- s ), where 

L 2 > s (R d ) = {u(x)\(l + \x\ 2 )^u(x) G L 2 } 
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and 

T-L m , s = {u{x)\D a u G L 2 ' s , < |a| < m}. 
Also, for any / G L 2 ' s , 

R ± (X)f = Ro(X)f - R^(X)VR ± (X)f, 

where Rq is the resolvent for the constant coefficient operator, P. As the notion of short 
range deals with compactness of the operator Z{u) = R(Vu), being short range requires 
sufficient decay assumptions at oo on V. Heuristically, it is required that the coefficients 
of V decrease as fast as an integrable function in \x\ and for each fixed xq, we have 

The reasons why these heuristics hold true are explored below, hence we forgo this 
analysis here and move on with the fact that V(x, D) is a short range perturbation as an 
assumption. Note that in the case explored below, V is Schwartz in x and is dominated 
by -P(£) as |£| — > oo. It is also important to note that while contour integration works out 
nicely in M 3 , the results presented here hold in any dimension where Rq and Rq are arrived 
at through a limiting procedure. 

The Agmon approach to the distorted Fourier transform is equivalent to the approach 
taken by the author. Namely, we define 

4±(x,£) = e<* - R*m[Ve ix< ](x). 

Then, the distorted Fourier transform is a map T± : L 2 — > L 2 such that 

(i) Ker(jF ± ) = L 2 d1 where L 2 d is the restriction of L 2 to the discrete spectrum of P. Similarly, 
we have L 2 , the restriction of L 2 to the continuous spectrum of P. Then, the restriction of 
T± is a unitary operator from L 2 onto L 2 , 

(ii) for any / G L 2 

CF±/)(fl = (2tt)-' lim / f(x)~M^l)dx in L% 

and 



(F ± f)(x) = (27r)-i lim f f(OMx,Od£ m L 

3^°° JKi 



where Kj is an increasing sequence of compact sets such that UjKj = IR d \ M for 
M(H) = {£ G R d \\£\ 2 is an eigenvalue for H} U 0, 

and 

(iii) If P c is the projection of L 2 onto L 2 , then 

(P c H)f = (KM P{ ^±)f 
for any / G D(H) where Mp^ denotes multiplication by P(£). 
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In addition, we have ||-P c /||l 2 = ||-7~±/IU 2 - in other words, we have a Plancherel theorem 
for our distorted Fourier basis. 

Now, |Ho2j . Chapter 14 arrives at the same conclusions using 

However, using the resolvent identity 

R{z) = R,{z){I + VR {z))-\ 

we will see that a formal iteration shows equivalence between these definitions for £ large. 
It is precisely this iteration we use below to get uniform bounds in £. 

6. Convolution Kernels 

In this section, we derive the integral kernel in M 3 for the inverse of the differential 
operator 

P, = -A-|£o| 2 

where we have set /i = |£o| for simplicity. This will be quite useful in deriving the distorted 
Fourier basis functions for more complicated operators belows. 

Specifically, given u, f : M 3 -> R, we find K^x, y) such that if 

i> = /, 

then 

u= K t ,(x,y)f(y)dy. 
To begin, we Fourier transform the equation to see 

(e - n 2 )u = f, 

hence 

u=F- i i(e-v 2 r 1 ]*f- 

So, 

K,{x ) y) = F- 1 [{e-^y 1 ]{x-y), 

if we can define 

G(x)=^- 1 [(£ 2 -/i 2 )- 1 ](x) 

in a meaningful sense. 



DISPERSIVE ESTIMATES FOR MATRIX HAMILTONIANS 



9 




53- Re 



Figure 2. The contour for computing the behavior of 



Without loss of generality, set // > 0. Initially, asssume that i^O, though this will be 
easily seen as a limiting case in the end. We have 



2 _ ,,- ' 

OO f2TT 





4:71 



,i\x\rcos(fi) 



o Jo (r-fi)(r + n) 
rsin(lxlr) 

-dr 



r 2 sm(9)d9d(j)dr 



\z\Jo (r-n)(r + fi) 

by first making a rotational change of variables where £3 — > , then using polar coordinates. 

Now, we are set up to use contour integration to find G(x). See Figure [2] for the contour 
over which we integrate. We call this contour T^. 

Then, we have from residue theory 



ze 



iz\x\ 



(z-fi)(z + fl) 



dz 



2m 



lie 



i\x\[i 



2/i 



me 



i\x\fi 



However, breaking T down, we also have 

7 iz\x\ pR 



ze 



(z-fi)(z + fi) 



, rsm(\x\r) , ixie 1 ^^ 

dz = 2i I . v ' ' ' dr H 

[r-fi)[r + fi) 2 



me 



-i\x\fi 
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Combining terms and taking R — > oo, we have 

4vr 2 cos(/x|a;|) 

G(x) = : : . 

\x\ 

This is valid for all x since the integral diverges as x — * 0. 

Using simple residue theory, taking the distributional conventions 

/(A) = /(A + iO) 



or 



result in 

(6.1) G ± {x 
To see this, define 



/(A) = /(A-<0) 
4 7r 2 e ±«MM 



x 



Now, we may make the same change of variables and do contour integration as above, 
though in this case we need not worry about avoiding ±|£ |. So, our contour F Ri0 is the 
hemisphere on the upper half plane formed by the real axis and the half circle of radius 
say R> //. The only residue in such a region would be given by z = \i + ie as z = /i — ie is 
outside r^o- For each e, we then have 

Att 

G+(x) = —j^e-^. 
\x\ 



Taking e — > gives formula (6.1) for G + . The analysis for G is similar. 

The above analysis is then easily seen to be equivalent to applying to the distributional 
connvention 

/(A) = ^ [/(A + iQ) + f(X - i0)] , 

namely the case where both residues lying on the real axis must be taken into account. 
However, since our eventual goal is to work with oscillatory integrals, for convenience and 
without loss of generality, we will work with the complex operator G(x) = G + (x). 



7. Distorted Fourier Basis 

Note that in the sequel, we take the convention that the soliton parameter is A 2 instead 
of A. This serves to remind the reader of the positivity of this parameter. The convention 
of A slightly simplifies the variational formulation, but has no impact on the linear analysis 
presented here. 
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We seek to understand the functions in the continuous spectrum of 7i by decomposing 
them using a distorted Fourier basis given by 

(7.1) (-A + A 2 - Vi)(-A + A 2 - V 2 )u io = (A 2 + \Zo\ 2 )\ , 

where u^ = e lx ^° + g^ and g^ is yet to be determined. 
From (7.1 ), 

[(-A + A 2 ) 2 -(A 2 + |£o| 2 ) 2 K = (-A + A a )V r 2 « & +y 1 (-A + A a -V r a )ti eb . 

Hence, 

(7-2) [(-A + A 2 ) 2 -(A 2 + |eo| 2 ) 2 ]^ = F (D (x)e 1 ^ + V(x,D)g eo , 

where 

V(x,D) = V 1 (-A + X 2 -V 2 ), 
and F£ (x) is a Schwartz function. Then, taking the Fourier Transform, we have 

m\ 2 + A 2 ) 2 - (l^ol 2 + A 2 )%, = F(£;£o) + (V>^ )(0, 

where 

fe?)(0 = \ 2 (v 2 + v 1 )*(g) + (\ev 2 )Hg) + (v 2 + v 1 )*^\ 2 g) 

+ m)*m-(V 1 V 2 )*(g). 



Given 



L, = [(ie| 2 + A 2 ) 2 -(|^ | 2 + A 2 ) 2 ] 

= [(iei + ieoi)(iei-ieoi)(iei 2 + 2A 2 + ie | 2 )], 



we have 



K±*F $0 + K±*(y{x,D)g (0 ) 



where 



and 



Lt = m\ + led ± - 1^1 t me + ^ 2 + i&m. 

Note that for simplicity we have omitted a small complex perturbation in the elliptic term 
(|£| 2 + 2 A 2 + |£o| 2 ) since it does not effect the analysis. 
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V 

Figure 3. The contour for computing the behavior of the fundamental so- 
lution in the limiting case. 



To explore further, we see in 



f del + u ± io)(iei - led t + 2\* + 
(lei + i&i ± *o)(iei - m t ami + KoDciei 2 + + ie | 2 



using the change of variables £1 



\x\ 



2ir i*7v 



JO 



R 



Then, we have 



3 rcos(0)|x| 



(r + |Co| ± z0)(r - |£„| T *0)(r 2 + 2A 2 + |Co 



-r sin(9)drd9d<f). 



Doing integration first in 9, then a contour integral, we have as in Section [6] that 



so 



7T 



l£o| 2 + A 2 



s ±*[a>[[€o| _ e -N\/iioP+2A5 



For simplicity, we take = i£^(:r) as the analysis for i£T will be similar. Then, we 

want to use an iterative argument to show that for mid to high range frequencies, these 
distorted Fourier bases exist in L 4 . It will become clear in the sequel why L 4 is chosen. 
Note that since near 0, K is bounded, we have K e L 3+s for any s > 0. In particular we 
show the following: 



Lemma 7.1. For the operator K defined in Equation (7.3), we have 
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Proof. We actually prove the result for 



K(x,y)f(y)dy. 



The proof for K will be essentially the same. 
Using distribution theory, we have for s G M 

Kl(x) = 5(x) 



4tt_ 2 
x\ 



As convolution operators, 



and 



K°:L 2 -^ L 2 (0(1)) 
K 2 : L 1 ^ L°° (OdW 1 )), 



hence we wish to wish to define K s in such a way as to preserve these estimates and such 
that k s is analytic for < Re(s) < 2 and continuous for < Re(s) < 2. However, after 
making a branch cut on the left half of the real axis, for s G R we have 



ll(iei + ieol)- ls (ie|-|eol)- l 7(OII^ -wju 



< 



and continuity on < Re(s) < 2 follows easily on a strip in the complex plane. For ana- 
lyticity inside the strip, it is clear any factors gained taking derivatives will be logarithmic 
and hence controlled by the polynomially decaying coefficients from Re(s). Hence, using 
complex interpolation 

K 1 : L A (Ofl&l"*)). 

□ 

For simplicity, we from now on write K instead of K^ Q . Now, we seek to analyze the 
equation 

(7-3) g i0 = Kf o * F & + Kf o * (V(x, D) 9(o ), 

In particular, we have the following: 
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Theorem 3. Let P(x, D) be a differential operator of the form 

P(x, D) = (-A + A 2 - Vi)(-A + A 2 - Va), 

where V\, V 2 G 5. Assuming that there are no eigenvalues embedded in the continuous 
spectrum [A 4 , 00), there exists g^ Q G L 4 such that Equation (7.3) is satisfied for u^ = 
e tx ^° + gf {x)- We have 

g^x) = K ± *[f (-^ ,\^\)}, 

where fo is smooth in x, £0, |£o|> & n d 

\{x) N dZfo\ < 1. 
Moreover, there exists a value M such that for £ > M , 







where 

(7-4) \(x) N 9&dSf{x,to)\ <\io\ 2 ~ lal , 

for any multi-indices a and (3, N > 0. 



Proof. The solution to (7.3) will be solved differently for large and small values of £o- hi 
particular, we use a Fredholm theory approach for the small frequencies and an iterative 
approach for the large frequencies. The analysis will be done using K + as the analysis for 
K~ will follow similarly. For simplicity, we set K = K + . 

To begin, let us take |£o| > M, where M will be determined in the exposition. Then, we 



.2 



solve Equation (7.3) using Picard iteration. For simplicity, let g^ = v. Setting v° = and 
Tu = V(K * u), we have 

v 1 = K{x)*[F io {x)e ix<0 } 

= K{x) * [(FUx)e lx <°) + (V(x,D)K(x) * (F fo (x)e^°))] 

= K{x) * [(F fo (x)e 1 **>) - (V, + V 2 )(X 2 + \Zo\ 2 )K(x) * (F io (x)e^) 

+ (V! + V 2 )K(x) * (F io (x)e ix <°) - (W 2 ■ VK{x) * {F io {x)e ix < )) 

- (Vi(x)\/ 2 (x) + AV 2 )K{x) * (F^x)e lx < )] 



v 



K(x) * [F^ (x)e ix<0 + V(x, D)v n ~ ! 

n-1 



A-(,-)=:= ^ V'"/-; : (.r) ( "- 

m=0 
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We wish to show that this iteration converges in L 4 . To see this, let u G L 4 . Note that 
\\K * V(x, D)u\\ L 4 < \\K * Vu\\ L 4 + \\VK * Vu\\ L 4 + \\AK * Vu\\ L *, 
where V, V, V G S. Then, 

\\K*V(x,D)u\\ L * < * ] , 2 A \\Vu\\ i + I^L-Hdyl- 1 ) * Vu\\ L * 

■So + A |4o| 2 £o + A 

so using the Hardy-Littlewood-Sobolev inequality and the bounds on K, we have 

\\K*V{x,D)u\\ L * < \£o\-HVu\\ L * 

^ l^or»||V r ||£a||«||r4, 

for some V G S. As a result, 

\\K*V(x,D)\\ L ^ L * <C\^, 
where C is determined by V\, Vi- If |£o| > C 2 , then 

\\K*V(x,D)\\ L *_ L < < 1, 
and the existence of g% G L 4 for 

(I-K*V(x,D))gt = gt 
follows from a contraction argument. In the notation from the theorem, we have C 2 = M. 

Now, for the smaller frequencies, we apply Fredholm theory. This approach also works 
for large |£o|> however the iterative approach gives us uniform bounds for all £o such that 
|£o| > M. Once differentiability in £o has been obtained, we will then have uniform bounds 
for all £o- However, we must be careful near £o = as K has a particularly challenging 
dependence upon |£ |- We explore this shortly, but first let us finish the existence argument 
for low frequencies. 

To begin, Equation (7.3) shows that 
(7.5) g io = K* (V(x,Z )e ix *°) + K* (V(x,D)g^, 

where 

V(x,D) = (-A + A 2 -V^Vi + V^-A + X 2 ) 

is a second order operator. 

Now, if K* (V(x, D)-) is a compact operator, we may use Fredholm Theory (see [Evans] . 
Appendix F) to say that either there is a unique solution to (7.5) or there exists a nontrivial 
u G L 4 such that 

(I - K *V)u = 0. 
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However, expanding the equation for u, we see this u is an embedded resonance and hence 



an embedded eigenvalue from |ESlj or Mar-spec] . As our spectral assumptions preclude 



the existence of embedded eigenvalues, the solution to (7.5) is unique. 
Let us now discuss the compactness. The operator itself is of the form 

r \ e i\x-y\M _ e -|*-2/lVK°l 2+2A2 l ~ 
2 r e *l«-v[[&l _ e -k-s/lv / l«"l 2+2A2 " 

7T — 



\x-ym \ 2 + X 2 ) 
x [(-A, + A 2 - V 1 (y))V2(y) + V&K-Ay + \ 2 )]v(y)dy. 

Hence, using integration by parts, we are concerned about the following two types of 
operators 



(1) P lU = j K(x-y)V(y)u(y)dy 

(2) P 2 u = ! K(x-y)V(y)u(y)dy, 



where V e <S. Of course, technically there will be terms with derivatives falling on K and 
V, however a brief calculation shows that these fall into the same class of operators as Pi. 
Indeed, by construction 

(-A-|£ | 2 )^ = 



and 

4vr_ 2 
\x\ 



;-A - %\ 2 )K = ^L[ e -Wv^i^] 



hence when all derivatives fall on K, simply by looking at —A— |£o| 2 + |£o| 2 we get reduction 
back to Pi or P 2 as K is a convolution kernel for an exact solution. 

We now need to prove 

Pi ■■ L 4 - L\ 

for i = 1,2. 

Assume that Uj — > w in L 4 . Since we are working in IR 3 , using duality and the properties 
of V, we have 

PiUj(x) — > as j — > oo 

for almost every x, where i = 1,2. By the uniform boundedness of weakly convergent 
sequences, the Hardy-Littlewood-Sobolev Inequality, and Holder we have, 

< a 
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for i — 1,2. Hence, there is a subsequence jk such that || Pji^JI^ converges. Therefore, it 
must converge to 0. As a result, the operator K*(V-) : L A — > L 4 is compact and there exists 

~ 4 4 

a unique ^ for all £ - Note that VK is compact from — > L5 using similar arguments. 

To discuss the continuous dependence upon £ 0! we need to study the functions g^ in 
more detail. In particular, we must have Vg^ smooth with repect to £o an d |£o- From the 
expression for g^ , we know that 

(I-K*(V(x : D).))gs = (I-P)g^ 

so 

g (o = (I-P)- 1 (K*(V(x,£ )e i *>)), 

where 

^(0 = [(-A-e 2 )(-A + 2A 2 + e 2 r 1 - 
From Fredholm Theory and the spectral assumptions, (/ — P)^ 1 is a resolvent which is 
uniquely defined. However, using the decay of V, we can write 

V = v ± v 2 , 

where |e c l x IVi| < 1, \e c ^V 2 f \ < \\f\\ W 2,oo given < c < c . The constant c is determined 
by the decay of V. Hence, using a resolvent identity, we have 

Vgt = Vi(J - V 2 KV X )- X V^K * (V(x,Co)e Ho )). 

Using the decay properties of Vi for % — 1,2 and the differentiability of if, for any £ we 
have V 2 KVi(z) is well-defined for z G C in a small neighborhood of |£ |- As a result, 

(I-V.KV^- 1 

is analytic with respect to z. Also, K is analytic with respect to |£| and £, t^e lx ^ is analytic 
with respect to £ and we see that (?£ depends smoothly on |£| and £. Using the resolvent 
identity 

f (x, = Ve ix < + V{1 - KV)- X K * {Ve lx <), 
the decay in x for f follows. 

For |£ | > M, let us return to the iteration scheme 

gl = K*\V(.,t )e*M], 

for n > 1. Assuming ^ = e tx ^°fo(x, £o, |£o|), we have 

fo = V(x,Z ) + e-^VK*(e^f ), 
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where by the mapping properties of K, choosing M large enough, this expression is valid 
in LJ for all |f | > M. 

We would like to better understand the regularity in x and £. To begin, let 

Then, we see 

(d x -iSo)u(x) = (^-^o)(^*[0(-,eo)e^ o) ])(x) 

+ J K{y)S x -rt»4> x {x-y,S Q )dy 
= J K( y )e^-^Mx-y,Zo)dy. 

From here, recognizing that e~ lx ^° cancels from 

and again using the mapping properties of K, we have 

ll^/oll i <c a , 

for all multi-indices a. Hence, fo G C£° fl Lf. Similarly, 

W(x) N dZM\ T i<CN„ 



for any N > using the decay in x of the operator V. 
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For the regularity in £, note that taking once again u = K * [e^'^°-*0(-, £o)]> we have 
(dt -ix)u = (^ o -za;)(^*[0(-^o)e l( -' 5o) ])(x) 

r-2 



+ % 



- ixj K{y)e^ x - y ^ct ) {x-y)dy + i / K{y)S v ^y<f>{y,^)dy 

+ f K{y)e^-y^(j>{x-y^)dy 
1 1 



K{y)e^y<f>{y,^)dy 



V^+2P [Col 

where we have used i£oe = d y e iy ^° and integrated by parts. As a result, 
for any multi-index /3, |/3| = 0, 1, 2, ... . Combining the above results, we have 

|^f/ (^OI<Ca,/3|e| 2Ha| , 

or /o G S* 2 , which gives (7.4). 

For the spatial regularity result, we once again use that the distorted Fourier basis 
satisfies the equation 

g i0 = K*(Fe ix ^) + K*(Vg i0 ). 

We have existence for g^ in L 4 , but we can take advantage of the structure of K * P in 
order to show improved regularity. Then, 

V<fe = (VK) * {Fe***°) + (VK) * (V^ ). 

Hence, we must explore the nature of (VK) * (V). Upon differentiating, we see 

(VK)(x-y) = 0(\x-y\~ 1 ), 

which means by a similar approach to Section [7], we get 

\\VgG>h*<c(\\F\\ L vt + 

To see this, we first use the Hardy-Littlewood-Sobolev inequality (see |Steinj ) with 7 = 1 



so 



1 _ 2 1 _ 11 
p ~ 3 + 4 ~ 12' 
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then Holders inequality such that 

\\Vg\\ L i2 <\\V\\ L *\\g\\ L *. 

Then, we can iterate this for all derivatives and using Sobolev embeddings, get continuity 
of all derivatives and hence smoothness. 

To prove existence for d^ g^ in Sobolev spaces, we must show that d^ g^ is defined and 
bounded in some space of functions. In this direction, we look at 

[(-A + 2A 2 + £ 2 )(-A-£o 2 )ko = F^ + Vg* 

and 

[(-A + 2A 2 + (£ + h 3 f){-A - (Co + M 2 )]^- = + Vgto+h,, 

where hj = hej and tj is the unit vector in the j-th coordinate. Hence, if we define 

then we must solve 

= O(h)(Fs + Fs + K*Vus )+V(v h ). 

We can write this as 

L(o[v h - 0(h)K * (K * (VgJ)] = 0(h)(G) + V[v h - 0(h)K * (K * (Vg^))}, 
where we have 

G = F (0 + F 6i -VK*(K*(yg 6i )). 
To see that G G L 4 , we need only see that 

\\VK*{K*{Vg^))\\ L A<oo 

since the other terms are dealt with above in the spatial regularity analysis. However, we 
have 

K * (K * •) : L 1 -> 

following analysis similar to the complex interpolation argument. Also, by moving all of 
the derivatives onto Pu, we see this is smooth. All we lack is nice decay, hence 

\\VK * (K * (Vgs ))\\ L * < \\K * (K * (V^ ))||l~ \\V\\ l *, 
for V E S as given in the description of P. From the Fredholm Theory, we know 




-O(l)K*(K*(y 9s0 ))\\ L *<C, 
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for C = C(£o)- However, given w G U L 4 a sufficiently decaying, smooth function, we 
have 

WvjWl* < C(l + \\wK*{K*{V gSo ))\\n) 
< C 

from Section [7], where C is independent of h. In this case, we have 

K*{K*{Vg (o ))eL°° 

using Holder's inequality, so we can take w = Thus, we can take the limit as h — > 

to see that derivatives in £ are bounded in weighted L 4 spaces. Iterating this process 
involves taking stonger weight functions at each step of the iteration. As a result, since 
V has exponentially decaying terms in x and Vg^ is well-defined in L 4 from the spatial 
regularity, we have the desired regularity in £ - 

Now that we have differentiability with respect to £oj 



% 0)j (j(-A + 2A 2 + £ 2 )(-A - £ 2 )k = Fe**> + Vg $0 
which implies 

L «o%o)^Co = d mj (Fe ixio ) + Pd (io)j g io 

- 2(6),(-A - £ 2 )^ - (&)(-A + 2A 2 + e 2 )^ - 

For higher derivatives in £ , we iterate this procedure. 



□ 



Remark 7.1. Note that the above analysis can also be done in the case where instead of 
L A we use L 2 ((x)~ s ) as in Agmon| . To see this, note that 

i < 



L 2 {{x)°)i 



where s > d, and 



l 2 {(x)-°) 



< 



where s > d. Then, we can go to the Sobolev norms to apply Hardy-Littlewood-Sobolev and 
use Holder's inequality in weighted spaces and the boundedness ofV± and V2 in weighted L 2 
spaces to complete the argument. 



Remark 7.2. As x — > 00, note that since V\, V2 G S, using Equation (|7.2|), we have 

■2 



71 



l£o| 2 + A 2 



\X\ 



which explains the choice of spaces L 2 ' s for x> \ in |Agmon . 
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8. Representation of the solution 



We present here a slightly different approach to the distorted Fourier transform, though 
the motivation comes from |Ho2j . 

Theorem 4. For V G S, there exists a distorted Fourier basis <p£ and correspondingly a 
distorted Fourier transform Q for the nons elf adjoint operator 7i, where 



Q±f 



4>f(x)f(x)dx. 



Similarly, there exists an inverse Fourier basis and correspondingly an inverse 

Fourier transform Q~ x for the nons elf adjoint operator Ti, where 



It follows that 



\Ql x \ 



< i, 

< i. 



These operators are not unitary, however 

||^± 1 ^||l2^ L 2 < 1 

and 

G±G±<I> = P c <\>. 

Before we prove the theorem, look at the operator 

ii 



L_L + 




for which we have the following self-adjoint realization 



LiL+Li 











Since 



LiL+Ll 



= (-A + A 2 - 7i)s(-A + A 2 - Vi - V 2 )(-A + A 2 - 
= (-A + A 2 - VJ 2 - (-A + A 2 - V l )W 2 {-A + A 2 - 
= L 2 _-LlV 2 Ll. 

This is a fourth order constant coefficient operator with a lower order perturbation. How- 
ever, the perturbation is no longer a differential operator. Ideally, by a similar analysis to 
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that in |Agmon], there exists a distorted Fourier basis, say such that 

LlL+Llut = (A 2 + £ 2 ) 2 tt£. 

i 

To prove this, we need to show Lt is a pseudodifferential operator of strong enough class, 
which we explore in the sequel. 

From Theorem pi we have = e %x( - + f$(x), = e lx( - + g^(x) such that 



H 2 



u 6 



u, 6 



where f$(x), g$(x) G L^., smooth in x and £, and 



7T 



ie i 2 +A 2 

as x — > oo. 

Formally, we would like to say 





i i 





LlL.Ll 



Llv e 



(A 2 + e 



2\2 



however as «£, t>g ^ L 2 , we must invest ige further. 

Before we begin, let us analyze the connection between and v^. For instance, 



L+(L_L + ui.) 
L_(L + L_t>g) 



Hence 



and 



In particular, we are interested in 

= (-A + A 2 -Vi)(e ia * + #) 
= (£ 2 + A 2 )e"'« + L_ 9i - V^, 

Then, C = (A 2 + £ 2 ), so 
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L -_\ = (\ 2 + er\, 

and 

A similar calculation holds for L + u^ = Cv%. 
Note also that if we look at the vector 




then we have 

i i 

To be more precise, we say that the operator LtL + L'i has a distorted Fourier basis given 
by Ug, then find an expression for the distorted Fourier transform of TCP C . This distorted 
Fourier transform will be defined via a distorted Fourier basis that will give the relationship 
between u^, and v%. The existence of must be proved since there is a lower order 
PDO perturbation instead of a differential operator. See [Ho2j. 

In order to prove L_ 2 is a PDO, we must use a result similar to one from |Ho4] . Chapter 
29. To this end, we refer to the following theorem given in |Ho4j : 

Theorem 5. Let X be a compact manifold, ^ a space of pseudo-differential operators and 
be the space of half- densities on X . Let P e fy™ hg (X; , Q2) be a positive, elliptic, sym- 
metric operator. Then, P defines a positive, self-adjoint operator V in L 2 (X, Q2. Ifm>0 
and a £ 1, then V a is also defined by a pseudodifferential operator in ^/^(X; , ), 
with principal and subprincipal symbols p a and ap a ~ 1 p s if p and p s are those for P. 

We seek to prove a slightly different version here: 

Theorem 6. Let P be a positive, symmetric, self-adjoint operator in \l/™^ 2 ' ) . Then, 
P defines a postive, self-adjoint operator V in L 2 (M. d ,M. d ) . If m > and a 6 1, then 
V a is also defined by a pseudodifferential operator in 2 ^(R d , WL d ), with principal and 
subprincipal symbols p a and ap a ~ 1 p s if p and p s are those for P. 

Note that since R G S, F(R) € S by the properties of the nonlinearity. Hence, we have 
the following: 

Lemma 8.1. The perturbation V± is short-range. 

We need to prove that given the operator, 

L_ = -A + A 2 - V x e S 2 , 
the new operator L a _ is a pseudodifferential operator for 06E. 
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Lemma 8.2. For an operator P, the resolvent R(z) = (P — z)^ 1 exists and is analytic 
for all z except the eigenvalues of P. Also, \\R(z)\\l2_^ L 2 is bounded by the inverse of the 
distance from z to the nearest eigenvalue. 



Proof. This follows from basic facts from spectral theory as discussed in |HSj . □ 
Theorem 7. The operator L a _ is pseudodifferential operator in the class S 2a for a G R. 

Before we prove the theorem, let us prove the following lemma from |Ho3j . 
Lemma 8.3. Let a G S m . If 

(8.1) \a(x,£)\>c\t\ m 
for |£| > C , then there exists b G S~ m such that 

(i) a(x,Ob(x,0-leS-\ 
(ii) a(x, D)b{x, D)-I G OpS' 00 , 

and 

[in) b(x, D)a(x, D)-I G OpS-°°. 

Proof of Lemma. First, let us prove that ( |8.1 ) implies (i). We can reduce this to the case 
where m = by looking at a(x,£)(l + |£| 2 ) _m/2 and b(x,£)(l + |£| 2 ) m/2 . 

Claim 8.4. If a 1; a 2 G S° and F G C°°(C 2 ), then F(a u a 2 ) G S°. 

Proof. Since the Reaj, Imaj G S* for j = 1,2, we may assume that ctj is real and 
F G C°°(R 2 ). Then, 



dF(a) = ydF_ d ^ 



dxj 2 -— ' dau 

J k 

8F(a) ^9F a 

-w - \^ a >" 

where d Xj dk G S°, dg.dk G S" 1 . Hence, it is clear the derivatives of F(a) decay as necessary 
for F(a) to be in S°. □ 

Hence, for m = 0, choose F G C°° so that F(z) = - for \z\ > c. Set b = F(a) G S° so 
a(x,£)b(x,£) = 1 for |£| > C. This proves {%). 

Using (z), we have that 

a{x, D)b{x, D) = I- r(x, D), r G S~\ 

Set 

6(x, £>)r(x, = h(x, D), b k G ^ m " fc , 
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so we can iterate out the error. Let b' be the asymptotic sum of the b^s, so 

a(x, D)b'(x, D)-I = a(x, D)(b , (x 1 D) - ^ b j( x , D )) ~ r ( x , D t e OpS~ k , 

j<k 

for every k. Then, we have (ii) replacing b with b' . Similarly, we can find a b" which satisfies 
(Hi). Note also that 

h' - b" = b'(I - ab') + (b"a - I)b', 
hence b' and b" are equivalent modulo S~°°. □ 

Proof of Theorem^ Since L_ is self-adjoint, we have that R(z) is defined and analytic for 
all z except at the eigenvalues of L_ . The L 2 norm of the resolvent can be estimated by the 
inverse of the distance to the set of eigenvalues. Now, since a < 0, we have by the spectral 
theorem 

/too 
z a R(z)udz, 
■ioo 

where the contour is slightly deformed near the origin to avoid z = and z a is analytic 
in the right half plane and equal to 1 when z = 1. Since = L^Lt, the distribution 
kernel of L°_ is an entire analytic function of a. 

To understand the behavior of the singularities, we construct a parametrix. Namely, 
since \L-(x, £)| > c|£| 2 for |£| > C, we have the existence of an inverse modulo S' -1 . Then, 
we can iterate that error, to find an inverse modulo S~°°. 

In particular, we have B z such that 

(P-z)B z = I-Q Zl 

where b z = F(P(x,£) — z), F(z) ~ 1/|^| for z large and Q z 6 Opi^S^ 1 ). Then, there is an 
E z given by the asymptotic sum 

oo 

J2B z (x,D)(Q z (x,D)y, 

j=0 

such that 

(P-z)E z = I-W Z} 
where W z G Op(S~°°). So, we have 

R(z) = E z + R(z)W z . 

Then, for a < 0, we have 

/ioo 
z a E z dz + T(a)u. 
■ioo 
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Here, T(a) should be analytic in a for a < 1. In particular, this remainder will be a well- 
behaved pseudo-differential operator using Beals' Theorem as discussed in |Beals] . From 
the composition of pseudodifferential operators, we have that 



Q z = J2 d t L -&0dZF(L_-z)/a\. 



a>0 

Hence, the terms of E z outside of compact set in phase space look like 
where q G S mk ~ K for some K > 0. 

_ i i 

Hence, there is a pseudodifferential operator representation of L_ 2 and thus L'i by 
multiplication by the operator. If p is the principal symbol of L_ , the principal symbol of 
L a _ will be F(p) where F(z) = z a for \z\ > C . 

□ 

i i 

Lemma 8.5. The pseudodifferential operator LJV\ + A + A 2 ) + LiV^Lf. is a short 
range perturbation. 



Proof. This proof should be similar to that in Lemma 8A The argument for the differential 

operator LJV\ + V\{— A + A 2 ) follows precisely as above. Hence, we focus only on the 

i i 

compactness and iteration arguments for the pseudodifferential operator, LlV^Lt- In what 

i i 

follows, let Tu = L^VL^K * (u). In particular, we need to prove: 

(8.2) LlVLle^ = e ix '^°V^ , for V^ G S, ||V£ ||x° ~ |£o| 2 > 

(8.3) ||^*T"(e^y eo )|U4 = O(\C \-^), 

(8.4) K*(LlVLb) : W 2A ^W 2 ' 4 . 



For (8.2), we have in the sense of distributions that 

Hence, since V G <S, 

L1VL1 = J P(x,0e l(x ' Xl ^V(x 1 ) J P(x 1 ,Ci)e ix ^5^ 1 )dCidx 1 d^ 

= J P(x^)e i{x - Xl ^V(x 1 ) J Pixu^e^dx^ 

= e ix t°V(x,£ ) + l.o.t., 

where V G S(x) and \V\ < £q precisely as in Section 7 This comes in particular from 

ii 

realizing that the principal symbol of L 2 _VL 2 _ is 

^ + X 2 -V l (x))V 2 (x). 



28 J.L. MARZUOLA 



The results (8.3) and (8.4) follow from the following theorem proved in [Stein] , Chapter 
VI. 

Theorem 8 (Stein). Suppose T a is a pseudo- differential operator whose symbol a belongs 
to S m . If m is an integer and k > m, then T a is a bounded mapping from W k,p — > W k ~ m,p , 
whenever 1 < p < oo. 

i ii 
Since Lt G S 1 and V G S°, we have LIVLt G S 2 , hence 

LlVLl : W 2 ' 4 -> L 4 . 
As V G S, we in fact have more than this. Define the symbol class 

st = { P \p g s m , \x a d?djp(x,o\ < Cajizr^}. 

In other words, we have the standard symbol class S m , where the symbol has rapid decay 
in x. Here, V G S®. Note that due to the properties of Schwarz class functions, we have 
for p G S mi and q G S m2 , 

pq, qp G 5 r mi+m2 

and 

qu : W m2 ' p -> L q , 

where 1 < p, q < oo. 



For (8.3), from the analysis in Theorem |3j we have 

(K* •) : W 2A . 



We have from (|8.2|) 



Then, 



K(x - y)LlV{y)Ll J K(y - z)^V io {z)dzdy\\ L , 
= \i Q \-l\\LlV{y)Ll j K(y-zy z ^V^)dzdy\\ Li 

using the fact that ^ 6 Sj and the mapping properties of K described in Theorem [3 
Iterating this procedure, we get the result. 
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1 1 
II T O T T T O _ II <S* 1 1 _ I. 
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|L»yL»u|| 4 < \\u\ 



By decay properties of V, we have 

i i 

\\XL1VL1U\\ L 4 X < \\U\\ W 2,4 + \\U\\ W 1A < \\U\\ W 2A. 

The inherent integration by parts is justified as V G S. Hence, by iterating this procedure 
and using properties of convolutions, 

[K*(LlVLl-)\ : W 2 ' 4 -> W 2 ' 4 ((-) N ), 
for any N eN. However, W 2,A ((-) N ) is compactly embedded in W 2,A , so (8.4) holds. 



□ 

i i 

Lemma 8.6. There exists a distorted Fourier basis, u^, for LtL + Lt with the aforemen- 
tioned smoothness properties. 



Proof. Apply the techniques from the proof of Theorem |3j applying (8.2), (8.3), and (8.4) 
when necessary. Once the compactness is established, the standard self-adjoint techniques 
are available to give 



(2tt 

Pr<t> 



■~ d 1 \F±<p\ 2 dx, 



where T± is the distorted Fourier transform associated to u^ Q and Po(£o) = (Co + * s the 
symbol for the leading order constant coefficient operator. 

□ 



Since 



J+ ^_ i 

LlL+Ll 





LJ 




' L 












o Ll _ 







L,L 



Ll 
LZ 



we have 



L 



Ll 



H 



Ll 
LJ 



Pcf = nw + ^Y\F±f, 



where J r ± is the distorted Fourier transform with respect to u^. Setting 



./ 



Ll 







Ll 
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we see 



LJ 


n 2 


Ll 


Ll _ 




LJ' _ 



pj = + yy\r± 



LJ 

o Ll 



Hence, 



H 2 (PJ) 



L_ 



rm 2 + \ 2 ?\(f± 



Ll 



or 



n 2 (pj)(x) 



f(Li^)(x)\(e + A 2 ) 2 | fZt(y)(LZ -*h){y)dydt 
[ j(L-J^)(x)\(e + X 2 ) 2 \ J ■ u^y){Llf 2 ){y)dydi. 
The inverse operations in these arguments are justified by the fact that 

L 2 = Ker(ft) Ker(^*) ± . 

We desire an oscillatory integral formulation for HP C - The continuous spectrum is 
spanned by the values ±(A 2 + £ 2 ) for all |£| > 0. Hence, we seek a diagonalization of 
the form 



HP C = Q- 1 



(\ 2 +e) o 
o -(A 2 +a 



Q. 



Using the above analysis for L'iL+Ll, we see that 

i(\ 2 + £ 2 )*FL~J (X 2 J^J^Ll 



Q = 



l 

71 

TTP, 



Q- 1 = 



1 

71 



-tLi^*(\ 2 +ey^ iLiF*{\ 2 +ey h > 



P- 1 F*T~\ 



where 



LJ 
Ll 



and 



T 



t(\ 2 +e^ (x 2 +ey^ 
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Note that we have for / = P c f 

H 

which is exactly what results from the decomposition. The resulting integral equation is 

-%J?_T*TJ?_f 2 



' fl ' 




iL-h 


. h _ 




_ -iL+fi _ 



HPJ 



So, since we have a pseudodifferential operator representation of Li, we could write HP C 
in terms of an oscillatory integral. 



Remark 8.1. We have now made precise the definition 
(8.5) fa = 



(8.6) 



z'(£ 2 + A 2 )Ll^ (^ 2 + A 2 )- 1 L!^ 

-tie + \ 2 )L~Jus (e + \ 2 y i Liut. 



where using the pseudo-differential analysis above, L_ 2 u^ is well-defined. 
Proof of Theorem \^ If 

" h 

h 



f 



e a ac {H), 



then 



Pf 



l\ 




" fl ' 


L~J _ 







e o- ac (H), 



where 



Ti 



LlL+Ll 






L^L+V 1 



Assume / G S, which we will relax later. Let ip be the PDO representation of P and 

$^(x) is the vector where both elements are the distorted Fourier basis function ^ for the 

i i 

self-adjoint operator L'iL + Li. Then, we have 

(0/XO = w,$ € > 
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where 

and $g is uniquely defined in the sense of distributions as 
where 

u = P(x, D)u^(x), 

and u (z S. Then, 

Similarly, we have 
where 

«£-! = p-^T- 1 )*^), 

where (T" 1 )* represents the adjoint of the multiplier matrix T _1 above. 

The modified Fourier transforms are in fact variations on the expansion involving the 
matrix Q. □ 

Note that since P(T~ X )* ,TP* G S°, the regularity properties of $£, $ _1 are the same as 
those of 0g as described in [3] with modifications to the explicit formulas. 

Note also the manipulations in the proof of [4] are valid in the sense of distributions, hence 
the assumption / e <S. However, as in the dispersive estimates below, similar estimates are 
seen to hold for less regular initial data through standard duality and limiting arguments. 

Corollary 8.7. As a result of the decomposition, we have a new proof of the fact that 

\\P c e m f\\ L i < 11/11*. 

Proof. This follows simply from mapping properties of pseudodifferential operators and the 
fact that the self-adjoint distorted Fourier transform is an L 2 isometry. □ 

Remark 8.2. Note that for convenience in terms of defining the resolvent, our result has 
been proved here only in IR 3 . However, using similar bounds developed in Ag mon] for higher 
dimensional resolvents, we expect that a result similar to that of\]\ holds in all dimensions 
and as a result similar estimates will follow below. The main difficulties presented would 
be a thorough discussion of the spectrum of TC as some of the known numerical techniques 
are unique to M. 3 . 



dispersive estimates for matrix hamiltonians 33 
9. Time Decay 

Using our distorted Fourier basis, we have that a solution to the problem 
(9.1) e iW P c = Q^e^Q^ 

for 

(a 2 + a o 



W 



-(A 2 + e 2 



The structure on Q allows us to do oscillatory integration in order to study the properties 
of e int . First of all, we prove Theorem [T] We will fix the notation K = K + . 

Proof o/[I} Using matrix notation, we have 

{Gim = J Mx)iJ( x )dx, 

where 



and (f)£ is given by (8.5). 

Looking at the integral representation, we have 

e itH Pj(x) = [ fc\x)e* w [ te(y)j(v)dyd£. 



Let x £ i ^ e a smooth, cut-off function chosen such that the iteration techniques in 
Theorem pi hold for £ G W 1 \ supp(x)- Then, take 

(9.2) e un Pj(x) = [ X mt^y m [ fav)j(v)dyd£ 



(9.3) 


+ 


J[i-x{mt^y tw 


/ <Pt(v)'$(y) d y d £- 

>y 


Hence, we must bound 












II = 


f[x(0 + 


(1-X(0)]^ 1 (^)e ± ^ 24 


Jy 


III = 


f[x(0 + 


(1 - x(OW x 'e ±u ^ 


\ gdv)^{y) d v d ii 


IV = 


f[x(0 + 


(l-xiOM 1 ^^ 


~ x2) [ gt{vMy)dyd£ 

Jy 



34 



J.L. MARZUOLA 



From henceforward, we work only with the term 



e i\x\\ 



\X\ 



from as the analysis for the exponentially decaying term will follow using simpler ver- 
sions of the methods for this case. Many of the techniques used are developed from the 



presentation in Schlagl . The challenge lies mostly in that <9?|£| is not bounded near 
for |a| > 2. Thus, we must be careful near the origin using stationary phase arguments 
since error terms require a minimum of two derivatives. A discussion of stationary phase 
complete with proofs is given in |EvZwj or [Stem]. Take the integral, 

where h(x) E C c °°, P{x) E C°°. Assume that d x P(0) = and d%P(0) ¥= °- Then ; the 
principle of stationary phase gives 



CLjT J , 



X ~ T~2 ^ , 

3=0 

where the asymptotic terms in the stationary phase expansion are given by 

dj = L j h(0), 

for V an order 2j differential operator as discussed in |EvZwj . 

Equation / is bounded using standard techniques of contour integration from the Linear 
Schrodinger equation. In particular, we have 



Before we investigate further, we recall some properties of the functions d^g^. From the 



expression (7.5) for g^ , we know that 

g io = K* (V(x,Z )e ix *°) + K* (V(x,D)g^ ), 

where 

Vg i0 = V X {I- ViKV^-^K * {V(x,£ )e* )). 

From Fredholm Theory and the spectral assumptions on Tl, (I — P^)^ 1 is well-defined, 
hence we can show that Vg^ is smooth in |£o| and £o- Also, K is smooth with respect to 
|£|, V^e %x( - is smooth with respect to £. As a result, as proved in Theorem pi g^ = K * f 



where f depends smoothly on |£| and £. Therefore, for £ near 0, we can take up to 3 
derivatives of the standard stationary phase operator 
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before we lose integrability in £. For £ large enough, from Theorem |3j we have 

9^o = K * f, 

where 

f = e**-t°f (x,Z Q ), 

where fo(x,£o) behaves like a symbol in S 2 . 

For ( |9.3 ), we use the principle of nonstationary phase and the principle of stationary 
phase in different regions. We have 

[1 - X(t)]%\x)e ttw J h(y)j(y)dydt, 

where 1 — \ is supported away from 0. In particular, we have integrals of the type 



i 



it 

where g and g are of the same form described above. Hence, we must bound the following 



II* = J[l-x(0M^y tie+X2) e-^(y)d^ 
III* = J[l-X®]e*e** +X %dvMvW> 

IV* = J[l-xmfa)e* e+X %(vMv)#- 

The bounds for e~ rf ^ 2+A2 - ) will follow through similar arguments. 
For integrals of type II*, we have oscillatory integrals of the form 

(q A) [ r i(-\z\\eo\-(x-z)to+tt*-y-to) fo( x - gjjgL f n 

{ } J N(A 2 +e 2 ) ^ 

Looking at the phase function, we have 

0(£o) = -\z\\Zo\-(?-z)Z +t£%-y-Zo, 
V^Mo) = 2ti Q -{x-z + y)-\z> h) 



vim = 



If we restrict £o to a region such that 



\z — y — x\ + \z\ 

u > V + 1, 
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then 0(£o) has no critical points. As a result, we can use the principle of non-stationary 
phase on this region with the decay properties of the function / to see we have decay like 
t~ N for any N. 

Let us hence assume that we are restricted a region 



I6>l< 



\z + y — x\ + \z\ 
2t 



+ 1, 



so cf) has at least one critical point. In fact, the critical point occurs where 

\z\ 



(9.5) 
For \z\ 



16 



z + y — x. 



y — x\ < 0, we have only 



(9-6) |£o| 
Otherwise, we have also 

(9-7) U 



\z — y — x\ + \z\ 
2t 



\z — y — x\ 
~2t 



As a result, all critical points occur on one of two spheres. Using (9.6) and (9.7), we 
have that if z, y and x are such that a critical point exists, that critical point is unique. 
Hence, we can define a cut-off function Xx, y ,z £ C^°(R d ) such that 



1 for |&| < I'-f-'l+H + M ? 
Ofor |,eo I > l'-g-*l+N + 



4 
M 
2 ' 



Let us assume that a critical point exists, say If |£q| < ^, the Hessian matrix is at 
least of rank 1 as £ <S> £ is a rank 1 matrix. So, there is at least one nondegenerate direction 
for £. After making an orthogonal change of coordinates bringing that nondegenerate 
direction to £i, using stationary phase on R, we have decay of the form 



||(9.4)| 



<t~K 



However, in the integral, we have rjr < r^n and |£q| > c > 0, so using the decay of fo in 
the overall decay is once again 

where the error is bounded by 



||(9.4)| 



sup I d? 

[a|<2 



fo(x - z,£p) 
(A 2 + £ 2 ) 



As fo & S 2 , this follows easily. 



DISPERSIVE ESTIMATES FOR MATRIX HAMILTONIANS 37 

For |£q| > the Hessian is nondegenerate. We can thus apply stationary phase in £ to 
get decay of the form 



IIHIU- <n 

where we have once again used the regularity of fo is x and £. Then, given the uniform 
decay of fo and boundedness in y and x, we have uniform boundedness with decay of type 
t~ 2 . The result for type III* follows similarly. 

The analysis for oscillatory integrals of type IV* is similar in that the phase function 
becomes 

0(£o) = Mlfol + (x-z)£ Q + t$ - \z Q \\£o\ - (y-zo)£ , 
V & 0(f o ) = 2t^ + (x-z)-(y-z ) + (\z\-\z \)-^, 

ISO I 

vim,) = 2u i + M^M(/,-^). 

l?o| l&r 

Hence, where critical points exist, we split up the regions of integration into |£o| > 
and |£o | < ^2*^°^ - Once again, we have stationary phase in full on the first region and 
stationary phase in at least one direction, coupled with the fact that A- < 2|Jbr Away from 
the critical points, we once again apply non-stationary phase. 

Let us now analyze (9.2). In particular, we have integrals of the type 

+ k(x)y tie+X2) I (e-** + 9c(y))iP(y)dydZ- 
; Jy 

Thus, we have to bound 



II** = j[ X (0]9^^)e U{e+x2) e-^, 
IIP* = Jm)]e lxi e^ 2+ ^g^y)d^ 



Once again, the bounds for e j *(£ 2+a2 ) will follow from similar techniques. 

For integrals of type II** and III**, we have an oscillatory integral of the form 

(9.8) fji-MM+tti-M M*-*) ^ 



38 J.L. MARZUOLA 

The phase function is 

0(£o) = -Ml£o|+^o-^o, 

V 5o <K£o) = 2tf -y- Ni^T, 

l?o| 



Isol |?o| 



Let us begin with an integral of type //**. After making the orthogonal change of coordi- 
nates £1 — > A and moving to polar coordinates in £, we need to bound 



/•27T p-K poo p —ir\zo\ 

< II / / / / x(r)^ T - r e i * r V ir|l ' |o ° sW 
JO •/ •/ o i \ z o I 

x /o(^ — ^o, r sin(#) cos(0), r sin(#) sin(0), r cos(#), r)r 2 sin(#)dr<i#<i0||i,°°. 



If we Taylor expand f in terms of (r sin(0) cos(0), r sin(6 l ) sin(0), r cos(#), r), then we can 
integrate in 0. In which case, all terms in the expansion with odd powers of cos(0) or sin(0) 
vanish under integrating out, leaving us with a function of the form 

/ (r 2 cos 2 (£),r). 

Integrating by parts in 9, we have 

p2n pn poo -ir\z \ 

< II / / / x {r)^^e ltr2 e- irlylcos(e) 
Jo Jo Jo \ z o\ 



x f (x - z ,rcos (6),r)r sm(6)drd6d(/)\\ L o 
< 



^ o X[ ) \zo\ \y\ 

x f Q (x - z ,r)rdr\\ L oo 

noo p -ir\zo\ nP ~ir\y\ cos(0) 

x dgfo(x — zo, r cos 2 (#), r)rdydzodrd9\\ 

'y Xl j M |y| 
x f (x - z ,r)rdr\\ L oo 

II f rvM^'y ^lglHg)) 
11 7o A xl j W \v\ 

x dgfo(x — Zo,r cos 2 (9), r)rdrd9\\ 
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since for n odd, 



1 • f 1 

e^ x x n dx = i J sm({ix)x n dx. 



'-l J -i 

Note that the boundedness in y and 9 is hence maintained after the integration by parts. 

Let us extend the region of integration in r to R. Due to the nature of the oscillatory 
functions involved, we experience no loss in doing so. Then, using the linear Schrodinger 
equation dispersion, we have 



I 7T**II < II / c itr 

\ 11 \\L^ ^ -II / 



■i 



oc 



X d r 



e -ir\z \ e ~ir\y\ _ e ir\y\ 



< 

~ , 3 



Foi \y\ 

i „ r° , f i 



x(r)fo(x - z ,r)r 



dr\\ 



t2 j-oo j \zo\\y\ 
--i 



x [F [x{r)fo(x - z ,r)r] (u + \z \ + \y\)(x - z ) 
- ^ [x(r)f (x - z ,r)r] (it + \z \ - \y\)(x - z )dy\du\\ L o C . 

From the estimate 

INU 1 ^ sup ||<9 Q u|| L i, 

\a\<d+l 

coupled with the facts that x £ Co°> fo £ C™, and / is rapidly decaying in x, we have 

||//**||l-<^||/||li. 

For the integrals of type III**, we immediately apply the linear Schrodinger estimate to 
get 

l|//nU~ < \ f\f x{O e ^^^h{x 1 ^Mf{y)dx 1 di\dy 
£2 j j \y — x\\ 

< 1 

r~j 3 5 

using once again the smoothness and decay of x, fo- 

The analysis for oscillatory integrals of type IV** is similar to that for type II**, except 
now we have no 6 dependence in the phase. Thus, we have phase functions of the form 

0(6) = -MI6l+*6 2 + MI6l, 



V^(6) = 2^o + (M-M)J 5 7 

l?o| 



\y 




x\ 


161 





V Co 0(6) = 2tl d + |fJ {I d - ). 
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At this point, it becomes convenient to move to polar coordinates in £. As a result, we 
have 

-27T n7V poo p -ir\z \ nP ir\zx\ 



\\IV**\\l~ < X(r. , , 

Jo Jo Jo \ z o\ \Zl\ 

x f (y — z , r sm(9) cos(0), r sin(0) sin(0), r cos(0), r) 
x / (z — 21, r sin(#) cos(0), r sin(0) sin(0), r cos(6 ) ), r) 
x r 2 sin(6)drd6d(f). 

Hence, we can first extend the interval of integration in r to R, then immediately integrate 
by parts in r to gain a factor of |. We once again apply the linear Schrodinger dispersive 
estimate to get 

t2 



Combining the above results, we have 
and 



IIOIU- <t-hy.<\\v 



d , 



IIQIUoc <r*\\<ii)\\ L x. 

Hence, the theorem follows. 

□ 

We now proceed to prove Theorem |2j 

Proof of\^ We proceed similarly to the proof of Theorem [TJ except now we must bound 
the following: 

i = |e- cN / x(0^ 1 (^)e ±i4(?2+A2) / MvMv)*u*t\> 

Jti " Jy 

ii = le-^i /[l-xlOl^W^^ f MvMv)*v*t\- 

Jti " Jy 

For II, we look at oscillatory integrals of the form 

/[I - x(£)](e*"* + 9s\x))e ±u{e+x2) [ (e** + ge(y))rKy)dyd£. 
Jti Jy 

Motivated by the principle of stationary phase in [EvZw] . define the operator 

r_ 

±2|£| 2 it' 

Considering the phase function as 0(£) = t£ 2 , it is clear 
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Then, let us take L M e 1 ^' in II and integrate by parts. Note, on the support of 1 — x(£)> 
£/|£| 2 is a bounded multiplier. A calculation shows 

|%d < I J ^e tlx - ym e^f (y,Ody\ 

vs + a 

P i[»-»[[€l _ P -|*-»lV^+ A2 

f - i/l 

+ I / j j e^/ (y,0^l 

< 

Using the regularity of f in y and £ an d continuing this calculation for d^g^, by applying 
the decay results from similar terms in [T] we see 



((x) + (y))\f (y,Z)\dy+ / \dMy,£)\dy. 



Now, for /, we need to bound 



bciZW* + 9e + Ke iy< + ge(y))d£. 

It is here our moments conditions become necessary. We wish to proceed similarly to case 
II, but now j|W is a singular multiplier. In fact, note that after integration by parts M 

times, the leading order operator will be on the order of |£| _2M . As a result, we arrive at 
the 2M moments conditions in (4.1). We have a gain in time decay using integration by 
parts in L, and since 

Ljg(0) = 0, 

for j = 1, . . . , 2M, there still no singularities near £ = where 

(9.9) £(0 = x{Z)fcH x ) I Mv)f(v)dy, 



and Lj is the order 2j differential operator resulting from the stationary phase-like argu- 
ments. Now, again we can apply the applicable results on oscillatory integrals of the terms 
77**, III** and IV** from the proof of Theorem [j] with new functions f^ 1 

f M (x - z,y,t, = \x\ A ' h y M2 m Ml ,M^ \m A ' h fo(x - z,t, 

defined on the support of x{0 where Mi + M 2 + M 3 = 2M. Using the moments conditions 
and the weighted integrability of /, the argument proceeds precisely as that near £ = for 
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the unweighted time decay case. Hence, under our assumptions we have 

l|e" cN / x(t)fc\*)e* w [ My)f(y)dyda^<t--^ M . 



□ 



Remark 9.1. In turn, (4.1) becomes our moments condition for the function space as 
defined by 

V2 = {0 G -P C 7Y 1 1 1 <^ 1 1 i^-A < 00, 111x1^011^2 < 00, condition \4-l\ is satisfied for j < A}, 



with norm 



v£ = {\m\ 2 H A + \\\x\ A (j)\\ 2 L 2y 



1 



These function spaces will be used in |Mar-nonlin] in order to find stable perturbations of 
minimal mass solitons. 



10. Dispersive Estimates 



From |Weinlj or Mar-spec| , we have H 1 = M <g> S where M is 2d + 4 dimensional set 
of functions that span the 4th order generalized null space at and S is the continuous 
spectrum. 

Since M is spanned by functions with exponential decay, we have for e M 



\e un <P\\ m <C(l + \t\ 



-c\x\ 



\4>(x)\dx, 



where c is determined by the exponential decay of all functions in M. 
Now, from [ESI] and [8] we have for 6 S, 

(10.1) ||e im 0|| L2 < C||0|| i2 . 



Lemma 10.1. Given Equation (10.1), we have 

||e* m 0||^ < CUWhi. 

Proof. For G S, we have 



\m < 

< 

< 
< 



m. 



\\He ltn <p\\ L 2 +C\\e 
\\e ltn H<f)\\ L 2 + C\\e itn <f) 
||H0|| i2 + C||e 4m 0|| i2 

\ H 2+CU\\ L 2 



< cm*. 

Hence, the result follows from interpolation. 



□ 
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In order to push through the contraction argument, we need various dispersive estimates 
from |BouWaj . We present the proofs here. 

Theorem 9 (Erdogan-Schlag,Bourgain). Let P c and P d be projections onto the continuous 
and discrete spectrum of TC respectively. Then, 

(i) \\e un P c <P\\ H1 < C\\(f>\\ w 

(«) „ _ 

(Hi) \\e it7i {P d <l>)\\ Hs < C(l + |t| 3 ) [ e- c W\4>(x)\dx 



) V m P4 




< 


V tH {Pc<P) 


Whs 


< 


\\e un (P d <P) 


Whs 


< 


\ a e un (P c 4>] 


lib 


< 


a e m (P d <P\ 


HI* 


< 



L* < C(l+|t| 3 ) J 



\H° 

\e~ clxl dx. 



Proof. Estimate (Hi) follows from the discrete spectral decomposition into a 4 dimensional 
generalized null space. The exponential decay is apparent from the properties of the eigen- 
functions. Estimate (v) follows similarly. 

For 4> G cr ac (7i), we have from Section M or |ES1] that 



For G c ac (H), we have 



||e im P c 0|| L2 < CU\\ L2 . 

? tH <PWm < Wne m <P\\ L 2+C\\e UH <P\\ L , 
< We Un n<f>\\ L 2 +C\\e itn 



L 2 



< \\H<l>\\v + C\\e ieH <l>\\ v > 

< ||0||^ + c\\<j>\\ L , 

< cu\\ H 2. 



This gives (i). A similar argument shows 

\\e un 0\\ H2s+ i < U\\ H ^ + \\e im ( l>\\ H ,s-t. 

Thus, by induction, we have (ii) for all positive integers s and hence by interpolation all 
s > 0. 

Let (f) G a ac (7i) and u = e ttH <j). Then, since 

iv t — Hv = 0, 
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then 

j t J \x\ 2al \v(x } t)\ 2 dx = 2Re(\x\ 2a v } v t ) 

= 2lm(\x\ 2a v,Hv} 
= 2lm(\x\ 2a v,Av) +0 (^J \v\ 2 e~ clxl 

< J \x\ 2a - 1 \v\\Vv\dx+ \\v\\ 2 2 . 
Using the following interpolation inequality 

II I rv I a 0) I 

II 2 



\x\ a -^\D^v\\\ L 2 < Wlx^vWl^ \\v\\l 



we have 



/N-Hiv^ < HN-HIUHIkl-IVHII. 



< \\\ x \ a \ v \\\ 2 L 2 a \\ v \\H«- 



Hence, using (ii) 



,! ''^n^iii^iiixrHii^-ii^iii+Hi,. 



dt 

Integrating, we have 

rt 

Ii-lm[o./]i ^ III-' I l"IN2 ••• / .III-' I <•«*'!!/- 





A a \v\\\l wm]) < \\\x\ a M\l+ [|IW B «WllS = II^III- + IK*)llr"]^ 



< IIN10ill2 + IIN a klll'lloo ([O)t]) / [UMc + U\\h)ds 

J 

< iiixii0iii^+£|iixiiHiii^ ([o> t])+^)(* 2a +t)iiixri0iii 2 2 . 

Hence, estimate [iv) follows. □ 

11. Strichartz Estimates 

From the above time decay, we can also prove the standard space-time Strichartz esti- 
mates for e lHt (f) where G P^H. We review the standard methods here as seen in |SulSul] . 
From henceforward, let us assume that we work on the subspace of functions contained in 
PcH. 

Theorem 10. For p and p' such that - + ^ = 1, with 2 < p < oo, and t ^ 0, the 

transformation e lHt maps continuously L p ' (R. d ) into L p (R. d ) and 

("•I) \\e mt <P\\L P <^r^U\\ LP '. 



2 p' 
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Proof. This result follows from the interpolation result presented in |BeLoj . 

Definition 11.1. The pair (q,r) of real numbers is called admissible if - — | — ^ with 
2 < r < when d > 2, or 2 < r < oo when d = 1 or d = 2. 



The following result proving Strichartz estimates is from Schlagl . 
Theorem 11 (Schlag). For every 6 L? and every admissible pair (q,r), the function 



JHt, 



belongs to L q (R, L r (R d )) nC(K, L (IR )), and there exists a constant C depending 
only on q such that 

(11-2) \\e mt <P\\ Lq(R , L r m) < C\\<P\\ L2 . 

Proof. Typically, one uses a duality argument when the operator e int is unitary. Namely, 

|(e^0,G) i2(Md+1) |< H\\ L 4G\\ Lq > L r>. 



To this end, write 



(e int 4>,G) L 2 m ds 



AH-s 



G(s)ds 



L 2 



< 



where 



AH-s 



G(s)ds 



L 2 



L 2 



JH- 



jn-s 



G{s)ds 



L 2 



, G{s)ds, 



-oo 
oo 



jn-t 



G(t)dt 



G(t)dt, 



AHt-s 



G(s)ds 



L 2 



< ll^ll^'Lr-' 



e m - s G{s)ds 



LiL r 



Using Equation 11.1 we have 



JHt-s 



G(s)ds 



/oo 
\\e int - s G(s)\\ Lr ds 
-oo 



< 



< 



1 



\t-s\ d ^-> 
1 



\\G(s)\\ Lr , ds 



\t — s\i 



\\G(s)\\ Lr , ds. 



Hence, using the Hardy-Littlewood-Sobolev Theorem with 7 



JHt-s 



G(s)ds 



< 



\\G\\ 
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However, for systems, this is not applicable. Hence, we must use the Christ-Kiselev 
Lemma |ChKi] . 

Lemma 11.2. Let X, Y be Banach spaces and let K(t,s) be the kernel of the operator 

K : L p ([0, T];X) — > L q ([0, T];Y). 
Denote by \\K\\ the operator norm of K. Define the lower diagonal operator 

K:LV([0,T];X)^LO([0,T\;Y) 

to be 

Kf(t)= [ K(t,s)f(s)ds. 
Jo 

Then, the operator K is bounded from L p ([0,T]; X) — ► L g ([0,T];Y) and it norm \\K\\ < 
c||K||, provided p < q. 

A perturbative approach originated by Kato is used. Define 



{SF)(t,x) 



Jo 



Then, 



ll^llipza < II^IU^i- 

Using the fractional integration argument from the unitary case, we have 

\\SF\\ L r L ><\\F\\ L , Lq j, 

where (r,p) is admissible. By Duhamel, we have 



e - itn P r = e- itHo P c - % 



e -i(t-s)H 0Ve -isHp ds _ 



Set V = MM~V, where 



Then, 



-i(t-s)H 



M 



Mg(s)ds 



(x)~ l - 
(x)- 1 - 



< 



e lsHo Mg(s) 



L 2 



where the last inequality follows from local smoothing. Applying the Christ-Kiselev lemma, 
for any Strichartz pair (r,p), we have 

,-t 



-i(t-s)H 



Mg(s)ds 
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Then, 



-un 



^/|L L? < 11/11^ + 



LlLl 



< 



L 2 ■ 



so we need 

M-^Ve- isH P c f 
Taking a Fourier transform in s gives 

/oo 
\\M- l V[P c (H - A - tO)P c ]- 1 PJ\\U\ < ||/||| 2 . 
■oo 

However, this follows from the smoothing estimate on Hq, plus the standard resolvent 
identity under the spectral assumptions on 7i. Hence, 

\\e- UH P c f\\ L r tLl < \\f\y. 

□ 
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